Domain Walls in Normal and Superconducting States of Iron Pnictides 
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The electronic and magnetic structures in the normal and superconducting states of iron pnictides 
are investigated by solving self-consistently the Bogoliubov-de Gennes equation. It is shown that 
strong electron correlations can induce domain walls, which separate regions with different spin 
density wave orders. At zero or low electron doping, 90° domain walls are formed while anti-phase 
domain walls are produced at higher electron doping. On the domain walls, larger electron densities 
are always present. The results agree qualitatively with recent observations of scanning tunneling 
microscopy and superconducting quantum interference device microscopy. 
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The discovery of a new family of layered superconduc- 
tors, i.e., the FeAs-based superconductors, could offer a 
new avenue to explore the mechanism of high temper- 
ature superconductivity [1-5]. Similar to the cuprates, 
the parent compounds of the FeAs-based superconduc- 
tors also possess the antiferromagntic ground state [4,5]. 
With increasing electron or hole doping, antiferromagn- 
tic order is suppressed and superconductivity appears in 
both the cuprates and the iron pnictides. However, dif- 
ferent from the cuprates, due to the the strong nesting 
effect between the hole Fermi surfaces around the T point 
and the electron Fermi surfaces around the M point , su- 
perconductivity and the spin density wave (SDW) orders 
can coexist in the electron-doped FeAs-based supercon- 
ductors [6,7]. Because each unit cell of the FeAs-based 
superconductors contains two inequivalent Fe ions, dif- 
ferent arrays of magnetic moments on Fe ions in both 
normal and superconducting states could lead to diverse 
magnetic structures and uncommon electronic properties 
[8,9]. 

Recently, in scanning tunneling microscopy (STM) ex- 
periments twin boundaries were observed in the nor- 
mal state of Ca(Fei_ x Co x )2As2 [10]. Across these twin 
boundaries, the a (b) axis of the crystal rotates through 
90°. This means that the modulation direction of SDW 
is rotated by 90°. In other words, 90° domain walls 
are formed at the twin boundaries. In Rcf. [11], su- 
perconducting quantum interference device microscopy 
(SQIDM) revealed that in the superconducting state of 
underdoped Ba(Fei_ x Co x )2As2 with x < 0.07, the dia- 
magnetic susceptibility is increased and the superfluid 
density is enhanced on the twin boundaries or 90° do- 
main walls. In another STM experiment [12], Li et al. 
also observed a 90° anti-phase domain wall in the parent 
compounds of iron pnictides, on which the local density 
of states (LDOS) is much higher than that in the interior 
of magnetic domains. Therefore, domain walls exist uni- 
versally in underdoped FeAs-based superconductors and 
affect strongly the electronic properties in the normal and 




FIG. 1: (Color online) Schematic of the tight-binding model 
Ho in Ref. [13]. ii is the nearest neighbor hopping between 
the same orbitals d xz or d yz on Fe ions, ti and tz are the 
next-nearest neighbor hoppings between the same orbitals 
mediated by the up and down As ions, respectively, and <4 
is the next nearest neighboring hopping between the different 
orbitals. 



superconducting states. 

In this work, we study the complex electronic and mag- 
netic structures in the underdoped FeAs-based super- 
conductors by solving self-consistently the Bogoliubov- 
de Gennes (BdG) equation. We start from the model 
Hamiltonian H = H + Hsc + H int . Here H is the two- 
orbital four-band tight-binding model proposed in Ref. 
[13], which describes correctly the characteristics of the 
energy band structure for the FeAs-based superconduc- 
tors [14-19]. The hopping parameters t\-t& in Ho are 
depicted in Fig. 1. 

The pairing Hamiltonian 



Hsc = X>i 
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where A^j,, is the pairing between the orbital /i (d xz or 
d yz ) on the site i and the orbital v (d xz or d yz ) on the site 
j, and cj is the creation operator of an electron with 
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spin a at the orbital fi on the site i. 

The interaction Hamiltonian Hi nt considered here only 
includes on-site Coulomb interaction U and Hund cou- 
pling Jtf. After taking the mean field treatment, H int 
can be expressed as [20] 

H int = U ^ ( n i^) n i^o- + {U -3 J H ) ^ ( n h^) n i"v 

+{U-2J H ) Y, < (2) 

where = 4^0-^ . 

We note that based on the model Hamiltonian H, 
the obtained LDOS [13,20], phase diagram [20], and spin 
susceptibility at different doping and temperature [21] 
for electron doped FeAs-based superconductors, and An- 
dreev bound states at negative energy inside the vortex 
core [22] for hole doped FeAs-based superconductors are 
all consistent with the STM [23-26] , nuclear magnetic res- 
onance [6,7], and neutron scattering experiments [27-30]. 

The eigenvalues and eigenfunctions of H can be ob- 
tained by solving self-consistently the BdG equation, i.e. 

El Hifj-jucr Ai/jjy \ / U- JVIJ \ _ / U iA1(J \ . . 

where Hi^ ua is the matrix element of H with spin a 
between the orbital /i on the site i and the orbital v on the 
site j. The superconducting pairing A^ji, = h(cintCji>± — 
Ci^c-jv^) in real space is associated with the eigenvalues 
E n and the eigenfunctions (u™ , v™ 9 ), and has the form 

n 

at temperature T . Here, fee is the Boltzmann's constant, 
and Vi^j,, is the pairing potential between the orbitals on 
the sites i and j . The corresponding local electron density 
reads 

* = EiK/JW + Km| 2 [1 - f(E n )}}, (5) 

where f{E n ) is the Fermi function, and the local mag- 
netic moment m, = \ Yj J&irf) ~ ( n if*l))- 

In order to interpret the complex domain wall struc- 
tures seen by STM experiments on the iron pnictides, we 
investigate the strong Coulomb correlation on Fe sites. 
In our calculations, we have employed the hopping pa- 
rameters in Ref. [13], i.e. t\ = 1, t 2 = 0.4, £3 = —2.0, 
and t\ = 0.04, and have chosen U — 4.8, Jh = 1-3, and 
Vi^u = 1.1 for fx = v and |i — j| = y/2, and for all other 
cases. Note that only the electron pairings between the 




FIG. 2: (Color online) The images of electron density m, 
magnetic order nti, and superconducting order A, at electron 
dopings x — 0.04 and 0.16. 

same orbitals on the next nearest neighboring Fe sites are 
considered. Such a choice of the pairing potential leads 
the superconducting order parameter to be s±-wave type. 

In Fig. 2, we present the zero temperature images 
of electron density rij, magnetic order n%i, and super- 
conducting order Aj = | J^j a ^i«jju a ^ electron dopings 
x = 0.04 and 0.16 on a 28 x 28 lattice with period bound- 
ary condition. From Fig. 2 (b) and (e), we can see 
that there exist obviously magnetic domain structures. 
Across the domain walls, the modulation direction of 
magnetic order rotates through 90° at x = 0.04 while 
the phase of magnetic order changes sign at x = 0.16. So 
the anti-phase domain walls, predicted previously in Ref. 
[8], are realized in the higher electron doped case. We 
observe that on both 90° domain walls and anti-phase 
domain walls, there are always higher electron densities 
m, as shown in Fig. 2 (a) and (d). Therefore, it is 
expected that superfluid density is enhanced on these 
domain walls, which coincides with the observations of 
SQIDM experiments [11]. However, the superconducting 
order parameter Aj has a larger magnitude on 90° do- 
main walls, but has a smaller magnitude on anti-phase 
domain walls (see. Fig. 2 (c) and (f)). 

In order to see clearly the variations of rij, mi, and 
Aj with distance from the domain walls, in Fig. 3, we 
give their values on the line y — 14. Fig. 3(a) and 
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FIG. 3: (Color online) Electron density m, magnetic order rm, 

and superconducting order A, on the line y = 14 at electron YIG. 4: (Color online) The images of LDOS p,( w ) with elec- 
dopmgs x = 0.04 and 0.16. tron dopingg x = QM and Q W at energies u = ± Q Q2 _ 



(d) show that more electrons are accumulated near all 
the domain walls. On 90° domain walls at x = 0.04 in 
Fig. 3(b) and (c), both mi and Aj have their maximum 
values. Oppositely, on anti-phase domain walls at x — 
0.16 in Fig. 3(e) and (f), almost vanishes, and A* has 
the minimum values. Therefore, except for the electron 
density, the magnetic and superconducting properties on 
the two kinds of domain walls are very different. 

We would like to mention that when x <^ 0.02, m and 
rrii have similar patterns with those at x = 0.04, except 
Aj = 0. Therefore, the 90° domain wall structure also 
exists in the normal state of the FeAs-based supercon- 
ductors, and is consistent with the observations of STM 
experiments [10]. However, we do not get the solution 
of complex 90° anti-phase domain wall seen in the par- 
ent compounds [12], which cannot be formed under the 
period boundary condition. 

Our calculations also show that with increasing elec- 
tron doping, the magnetic order m, is gradually sup- 
pressed and finally vanishes at x ~ 0.2, which is larger 
than the experimental value. The main difference be- 
tween the theoretical results and experimental data could 
be due to the fact that a strong Coulomb interaction 
U leads to renormalization of the hole Fermi surfaces 
around the T point and the electron Fermi surfaces 
around the M point in the FeAs-based superconduc- 
tors, which enhances the nesting effect between the hole 
Fermi surfaces and the electron Fermi surfaces. How- 
ever, by adjusting suitably the hopping parameters t% 
and £3 in Hq, which determine the sizes and shapes of 
the Fermi surfaces, the nesting effect could be dimin- 
ished and the experimental value can be obtained. When 



x <^ 0.08, the 90° domain walls always exist while the 
anti-phase magnetic domains show up at x >^ 0.15. 
When ~ 0.08 < x <~ 0.15, SDW and superconductiv- 
ity uniformly coexist. We also note that the solution of 
anti-phase domain wall structure is a metstable state in 
the above range, which has a slightly higher energy than 
the ground state. 

Now we calculate the LDOS on and near the domain 
walls in order to compare with the STM experiments. 
The expression of LDOS at energy u on the site i is 

= ^E{l<Mtkl 2 ^-o;) + |^ k | 2 ^nk + a;)}, 

n/j.k 

(6) 

where N is the number of wave vectors k, (wJutki u ii4k) 
and -EVik are the eigenfunctions and eigenvalues of the 
Fourier transformed BdG equation, respectively. In cal- 
culating pi(id), we have taken the delta function 6(x) = 
T/ir(x 2 +r 2 ) with quasipartical damping T = 0.005, and 
a 30 x 30 supercell is used. 

Fig. 4 shows the LDOS images on a 28 x 28 lattice 
at different energies and electron dopings. It is obvious 
that when uj = ±0.02, pi(u>) has the maximum value on 
the 90° domain walls at x — 0.04, but has the minimum 
value on the ant-phase domain walls at x — 0.16. Here 
we also omit the LDOS image with x = 0.0, which is 
similar to that of x = 0.04. 

However, the LDOS images change with energy uj. In 
Fig. 5, we give the energy dependence of the LDOS 
Pi(u>) at the sites on and near domain walls with dif- 
ferent electron dopings. In both x — 0.0 and 0.04, obvi- 
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FIG. 5: (Color online) Energy dependence of LDOS Pi(oi) 
with electron dopings x = 0.0, 0.04, and 0.16 at the sites on 
and near domain walls. The arrows point to the coherence 
peaks with maximum superconducting order parameter A, in 
the panels of x = 0.04 and 0.16 



ously, pi(w) at (14, 14) on the 90° domain wall is always 
larger than that on the other sites when ui € (—0.18, 0.1). 
In contrast, pi(uj) at (1,14) on the anti-phase domain 
wall is always smaller than that on the other sites when 
uj > -0.13. We note that when u £ (-0.2,0.2), the 
curves of pi(io) with x = 0.0 resemble those measured 
by the STM experiments in the parent compounds [25]. 
The coherence peak at positive energy is higher at both 
x = 0.04 and 0.16 due to the coexistence of SDW and 
superconductivity [13,20]. The asymmetry of the coher- 
ence peaks was also observed by the STM experiments 
[23-25]. 

In summary, we have studied the electronic and mag- 
netic properties in the electron underdoped iron pnic- 
tides. Due to strong electron correlations, the domain 
walls are formed at the twin boundaries producing be- 
low the structural transition of the parent compounds. 
The existence of the domain wall structures leads to the 
nonuniformity of the electron density and the supercon- 
ducting order parameter in real space. Therefore, in- 
homogeneity of superconductivity is intrinsic in the un- 
derdoped iron pnictides. The 90° domain walls and the 
supercurrent properties on them have been confirmed by 
the STM and SQIDM experiments. However, the anti- 
phase domain walls at the higher electron doping are not 
experimentally reported yet. We hope that such a mag- 
netic structure could be verified by future STM experi- 
ments. 
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